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In this chapter, we show two fascinating topics lying between quantum infor- 
mation processing and statistical mechanics. First, we introduce an elaborated 
technique, the surface code, to prepare the particular quantum state with ro- 
bustness against decoherence. Interestingly, the theoretical limitation of the 
surface code, accuracy threshold, to restore the quantum state has a close con- 
nection with the problem on the phase transition in a special model known as 
spin glasses, which is one of the most active researches in statistical mechanics. 
The phase transition in spin glasses is an intractable problem, since we must 
strive many-body system with complicated interactions with change of their 
signs depending on the distance between spins. Fortunately, recent progress in 
spin-glass theory enables us to predict the precise location of the critical point, 
at which the phase transition occurs. It means that statistical mechanics is 
available for revealing one of the most interesting parts in quantum informa- 
tion processing. We show how to import the special tool in statistical mechanics 
into the problem on the accuracy threshold in quantum computation. 

Second, we show another interesting technique to employ quantum nature, 
quantum annealing. The purpose of quantum annealing is to search for the 
most favored solution of a multivariable function, namely optimization prob- 
lem. The most typical instance is the traveling salesman problem to find the 
minimum tour while visiting all the cities. In quantum annealing, we introduce 
quantum fluctuation to drive a particular system with the artificial Hamilto- 
nian, in which the ground state represents the optimal solution of the specific 
problem we desire to solve. Induction of the quantum fluctuation gives rise to 
the quantum tunneling effect, which allows nontrivial hopping from state to 
state. We then sketch a strategy to control the quantum fluctuation efficiently 
reaching the ground state. Such a generic framework is called quantum an- 
nealing. The most typical instance is quantum adiabatic computation based 
on the adiabatic theorem. The quantum adiabatic computation as discussed in 
the other chapter, unfortunately, has a crucial bottleneck for a part of the op- 
timization problems. We here introduce several recent trials to overcome such 
a weakpoint by use of developments in statistical mechanics. 



March 2, 2013 16:21 



WSPC - Proceedings Trim Size: 9in x 6in O2'0hzeki 



2 

Through both of the topics, we would shed light on the birth of the inter- 
disciplinary field between quantum mechanics and statistical mechanics. 

Keywords: spin glass; phase transition; topological error correcting code; quan- 
tum annealing 

1. Introduction: Statistical Mechanics and Quantum 
Mechanics 

Quantum mechanics is a method partially based on the concept of the 
probability in outputs of the measurements on the physical state. This 
is because the physical state is allowed to be expressed by superposition 
of the possible eigenstates with the probability amplitude, which can be 
determined by the Schrodingcr equation. Readers might feel uneasy since 
such a probabilistic phenomena seems to be problematic and uncontrollable. 
However, as shown in this chapter, several techniques have been designed 
and proposed for realization of the stability in quantum systems. 

The quantum state is fragile, sensitive to noisy environment effects 
and continues to change. This peculiar destruction of the quantum state 
is known as decoherence. If we control to maintain the original state, we 
have to cure the quantum state suffering from undesired errors due to de- 
coherence. The technique to remove these errors is the quantum error cor- 
rection, which occupies the first part of this chapter. To fix successfully 
the quantum state, we need to prepare the redundant degrees of freedom 
to restore the original property and check where and how errors exist. We 
must thus design and deal with a many-body quantum system, which is 
known as an intractable problem. However it is worthwhile to strive such a 
difficult problem. Once you obtain an ingenious way to keep the quantum 
state, you find a pavement toward the realization of the quantum informa- 
tion processing, since you can encode simultaneously multiple information 
in the single quantum state by use of superposition. We here introduce a 
quantum error correcting technique, surface code, 1 which skillfully use the 
concept of the topology to encode the original state in quantum many- 
body system. The surface code has a remarkable feature closely related 
with the main concept of this chapter. The accuracy threshold, which rep- 
resents the theoretical limitation to restore the original quantum state by 
the surface code, is related with the phase transition in the special model 
in spin glasses. 2 Spin glass is a disordered magnetic material, which shows 
a peculiar behavior with extraordinary slow relaxation toward equilibrium 
state in a low temperature. 3-7 The complicated interactions in spin glasses 
spoil several methods used for systematic analyses in statistical mechan- 
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ics. The identification of the critical point in spin glasses is used to be an 
intractable problem for long days. Fortunately, the recent development of 
theory in spin glasses yields a systematic analysis for its precise location. 
It means that the specialized tool to analyses in spin glasses is available for 
the problem in quantum computation. We show the fascinating connection 
between two unrelated topics, while reviewing the systematic approach to 
derive the accuracy threshold. 

By using such a technique to maintain the quantum state as introduced 
above, we can prepare a coherent state. Then, how do we use superposition 
of the quantum system? Superposition enables us to evolve the quantum 
system in a parallel way over the various possible states. A famous algo- 
rithm given by Shor 8 is also based on the simultaneous search over all the 
possible candidates. In the second part of the present chapter, we introduce 
a generic technique by use of quantum nature in this chapter, called quan- 
tum annealing. 9-16 Quantum annealing is an active use of superposition in 
order to obtain the most important result by searching all the candidates. 
Its attractive feature is simplicity and generality in applications. This is 
useful to solve several particular issues known as the optimization prob- 
lems, in which we desire to find a minimizer or maximizer of the given 
multivariable function. 17,18 The Shor's algorithm uses a skillful technique 
in order to efficiently provide a desired answer. On the other hand, quan- 
tum annealing basically takes a simpler way only by controlling quantum 
fluctuations. The most typical procedure of quantum annealing, quantum 
adiabatic computation, is superior to its classical counterpart, simulated an- 
nealing, 19 ' 20 owing to usage of quantum fluctuation. However a bottleneck 
of quantum adiabatic computation is also revealed in the application to the 
particular optimization problems. This weak point can be understood thor- 
ough a special mapping from the classical stochastic time evolution to the 
particular quantum dynamics. We are again between statistical mechanics 
and quantum information processing. Therefore, in this chapter, we will go 
further to overcome the bottleneck of quantum adiabatic computation by 
use of this fascinating relationship. We show several attempts for surmount- 
ing the obstacle involved in quantum adiabatic computation by employing 
several alternative strategies from statistical mechanics. 

Statistical mechanics uses the probability similarly to quantum mechan- 
ics, but is capable to predict a definite behavior in the future in a large- 
number of components called as thermodynamic limit. The mathematical 
background of statistical mechanics, large deviation, provides the ability to 
give a definite answer even by use of the probability. Observant readers can 
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find out the key point in this chapter. That is to deal with a large num- 
ber of components. Turning on our eyes on substances around our daily 
life, they shows stable appearance. Material consists of many components, 
atoms, molecules and their mixtures, which should follow the rules of the 
quantum mechanics. The stability of these macroscopic systems comes from 
the particular property of the collection of many components. However, in 
many body systems, once if you tune the external parameters as temper- 
ature, pressure, and some fields, they can eventually change their uniform 
appearances. For instance, increase of temperature makes change of solid 
into liquid and gas. This phenomena is known as the phase transition. 
The phase transition involves singularities in the behavior of the physical 
quantities. The both of quantum error correction and quantum annealing 
shown in this chapter, suffer from this peculiar behavior. The failure to keep 
the quantum state by the surface code and the decay of the performance 
in quantum adiabatic computation come from the properties of the phase 
transition. Therefore, by dealing with the problems on the phase transi- 
tion, we study the fascinating interdisciplinary connection between statical 
mechanics and quantum information processing. 

Recent developments in manufacturing and manipulation of the quan- 
tum mechanical system enables us to prepare a large number of components. 
In this sense, the topics we deal with in this chapter must be valuable for 
developments in actual applications in the future. 



2. Training: Statistical Mechanics 

For unfamiliar readers with statistical mechanics and many students to 
understand the essential parts in this chapter, let us get back to the starting 
point of statistical mechanics. 

At first let us recall the concept of the probability. What was the prob- 
ability? If you have a non-tricky coin, you can say that the probability you 
can find a face (head or tail) should be a half. What does it mean? 



2.1. Student's misreading point: probability is... 

We use an artificial variable with a binary S = ±1 to represent the coin 
state as head (+1) or tail (—1). A non-tricky coin takes the probability as 
Pi (S = 1) = 1/2 and Pt(S = -1) = 1/2. We also can deal with a tricky 
coin with a biased probability as 

exp^S) 
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where K is the strength of bias, and Z\{K) is a normalized constant explic- 
itly given as 2 cosh K . By use of the above probability, we can evaluate two 
characteristic quantities of the probabilistic system. The first quantity is 
expectation expressed as m, which indicates tendency of the coin whether 
head or tail, 

m= SP 1 (S)=timh(K). (2) 
s=±i 

It readily shows that the coin state tends to be head if we increase the bias 
of the tricky coin K — > oo. Another is the variance defined through the 
square of the difference from the expectation as 

a 2 =^(5-m) 2 Pi(5) = l-tanh 2 ^). (3) 

s 

If we increase the bias K, the variance can vanish. However this is the case 
you know that this coin must be head without any probabilistic factors 
because P^S = -I) 0. 

Then, if you once glance at the result of the single trial, can you judge 
whether its behavior is closely related with the obtained value of m or not? 
The answer is "No". The value of m only gives "tendency". The actual 
results will fluctuate around the expectation in repeating observations. We 
must repeat the coin games many times. 

The above is the usual statement around the explanation of the proba- 
bility. Since students, who studied the probability theory or other related 
topics, do not find out its genuine meaning, they often say that "proba- 
bility" is obscure, abstract, and difficult!!, and conclude that it is not a 
well- understandable concept. Unfortunately, quantum mechanics and sta- 
tistical mechanics with the concept of probability might be considered as 
non-reliable and non-deterministic to predict the future-coming behavior of 
the system we deal with. I would like to ask you why the probability for 
a non-tricky coin was a half. The students must recognize what the expec- 
tation express in our future. As you know, the expectation is not relevant 
for understanding the output in a single coin game. Let us play the game 
repeatedly. Then... 

2.2. Probability describes... a certain behavior 

In order to clarify the concept of the probability, in particular expecta- 
tion, let us consider to accumulate the results of AT-time observations after 
flipping coins. In all the trials, their results are characterized by the individ- 
ually and independent distribution, which is same as in Eq. (I). Gradually 
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approaching the specific example of use of probability in statistical mechan- 
ics, we take a simple model of the magnetic material, which is the typical 
example of studies in this field. 

Let us introduce the Ising variable, which describes the magnetic mo- 
mentum of spins in the magnetic material, often simply termed as "spin" 
or "Ising spin". The Ising spin takes only two integers as ±1. We here 
regard the behavior of TV-independent spins as the nipping coins in the N- 
sequential times. In statistical mechanics, we choose the specific probability 
distribution depending on the conditions of the system under consideration 
in order to calculate the expectation. Notice that, in statistical mechanics, 
the probability distribution can be given " a priori" similarly to the above 
case of a non-tricky coin. In the previous simple case, we could choose a 
half as the probability for the coin since it was reasonable. Generally we do 
not know the precise structure of the probability distribution for various 
events. On the other hand, statistical mechanics provides us with the proba- 
bility distribution of the behavior in equilibrium state. For instance, for the 
equilibrium state without change of the number of components, we use the 
canonical distribution characterized by the Hamiltonian and temperature. 

P(S 1 ,S 2 , ■■■,S N ) = —L- cxp (-/3H(S U S N )) , (4) 

Zn{P) 

where we define the inverse temperature (3 — l/fc^T (often ks = 1) and 
Zn(P) is called as the partition function. The Hamiltonian describes the 
energy depending on the microscopic state of the system. The temperature 
controls thermal fluctuation to drive microscopic degrees of freedom in the 
system. If (3 — > 0, the spins are not stable due to strong effect by thermal 
fluctuation. On the other hand, in the low temperature (3 — > oo, the system 
is settled into a lower energy state. The single spin favors parallel direction 
to the magnetic field since its energy can be described by the Hamiltonian 
H(S) — —hS, where h is the strength of the magnetic field. Then the 
joint probability of the N spins (TV-time coin games) can be written by the 
canonical distribution as 

1 N 

P N {S U S 2 ,--- ,S N ) = — — TT exp(fihSi), (5) 
Zn((3) f \ 

where 

Z N ((3) = (2cosh((3h)f . (6) 

The partition function plays a roll of the characteristic function. Its loga- 
rithmic function is called the free energy —(3Fn{(3) = log Zn(/3). We often 
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use the thermal average by use of the canonical distribution in statistical 
mechanics written as 

(0) = J20({Si})P N (S,S2,---,S N ). (7) 
{s.} 

where O denotes the physical observable. 

Let me ask the following question. If you have huge number of the 
outputs in an experiment, how do you interpret the result. Usually we use 
the average over all the realizations to simplify the outputs. This simplified 
picture is often called as coarse graining in statistical mechanics. What we 
must take care of is that the average (sample mean) over all the realizations 

1 N 

and its variance 

1 N 

VN = ^Y,^- m N) 2 - (9) 

i=l 

Statistical mechanics gives a definite future of the above average for the 
many-body system in terms of the expectation and variance estimated from 
the specific distribution function. 

Below let us sec the connection between quantities given by the empirical 
measure and those by the a priori distribution function. 



2.3. Large deviation property 

We do not care the detail on the system with the large-number components. 
For instance, for the above N spin systems, it is enough to characterize the 
magnetic property by the average of the N Ising spins. This is called as 
the magnetization, which indicates the magnetic strength of the material. 
Statistical mechanics predicts such significant quantities in the macroscopic 
scale by averages over all the components, coarse-grained quantities. If we 
consider the large- number limit N — > oo, the averaged values can corre- 
spond to the expectations estimated by the probability distribution under 
a certain condition shown below. It is convenient to change the expression 
of the joint probability (5) into the form of the probability of the average 
as 

P N (m) = S [ m ~ jj 2 S i J Pn(Si,S 2 , ■■■ , S N ). (10) 

{S,} V i=l / 
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By use of the integral form of the delta function as 



/ioo 
diexp(xx), (11) 
-ioo 

we can obtain 

P N {m) oc f dmcxp<Nmm+ (-m^Si + l3h^Si\\ . (12) 
{Si}"* I \ i=l i=i J) 

Before performing the integration, we sum over the spin variables. 

P N (m) oc J dm exp{Nmm + N log (2 cosh - m))} . (13) 

The integrand is proportional to N. In that case, we can apply the saddle- 
point method to the above integration. When we consider the infinite limit 
of N, the integral is given by the maximizer to* of the integrand as 

P N (m) = exp {-Nf(m, m*)} . (14) 

Here we define the pseudo free energy /(to, to) as 

/(to, to) = —mm — log (2 cosh (f3h — to)) . (15) 

The saddle-point equation for to, which gives the maximum of the integrand 
(or the minimum of the pseudo free energy) , is 

df 

-4r = -)• to = tanh (fih - to) . (16) 
am 

The maximizer of the integrand is given by to* = fih — tanh _1 (TO,). By use 
of this maximizer, we obtain the probability of the magnetization as 

Pat(to) = exp | TV (^rnPh - TOtanfT 1 to + log — 2 | . (17) 

In Fig. 1, we describe the behavior of /(to, to*). The maximum is located at 
to* = tanh(/3/i), which can be verified by derivative with respect to to. If N 
increase, the probability away from to* = tanh(/3ft.) decrease exponentially. 
The maximum indicates the value of the expectation since 

Moreover, if we evaluate the variance, we reach 

- (j, f>) = V-W 1 = ^ " tanh2 ^) • ( 19 > 
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Fig. 1. Behavior of f(m,m*) with J3h = 1.0 (solid), (thick symmetric dashed 
curve), and —1.0 (dotted curve). 

Therefore the magnetization can take a "definite" value m* = tanh(/3/i) if 
N — > oo. This property is known as the large deviation principle in math- 
ematics. In the above simple case, we treat the individually-independent 
distribution. Fortunately, it is confirmed that the large deviation princi- 
ple holds in relatively wide varieties rather than individually-independent 
distribution. Statistical mechanics plentifully uses explicitly and implicitly 
this fruitful property. It means that, in many-body system, we can predict a 
definite future even by use of the probability as classical mechanics with the 
deterministic equation. In other words, statistical mechanics predicts the 
average value for the physical observables of the system with large number 
components from the calculation of the expectation by use of an a priori 
distribution function such that the canonical, micro canonical, and ground 
canonical ensemble. 

2.4. Mean-field analysis 

The probability can give a definite future expressed by the expectation 
value, if we are interested in many-body systems by the support of the 
large deviation principle. However many degrees of freedom in material 
strongly interact with each other. The spin is not an exception. Differently 
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from the above simplified case, we must take care of the effect through 
magnetic-dipole interactions to more precisely understand the behavior of 
magnetic material. We exemplify its simplest model to deal with effects of 
interactions known as the Ising model, whose Hamiltonian is 

H = -J^SiSj. (20) 

(ij) 

The subscript stands for the index of the site, where the spin is located. 
The summation is taken over the nearest- neighboring pairs of spins. The 
geometric property as the locations of the spins is closely related with the 
structure of the magnetic material. As precisely as possible we analyze an 
actual behavior of magnetic material, we must consider a finite-dimensional 
structure of atoms. However it is difficult to straightforwardly perform non- 
trivial analyses on the finite-dimensional Ising model. We take several ap- 
proximations or numerical simulations for obtaining meaningful results of 
the finite-dimensional Ising models. The exactly solvable examples without 
any approximations are found in one and two dimensions but rare excep- 
tions. 

The simple but qualitative approximation often taken as the first trial 
is the mean-field analysis. It enables us to obtain a meaningful picture 
to describe a peculiar behavior of the many-body system. The probabil- 
ity distribution of the degrees of freedom should be correlated with each 
other in the case of the existence of the interactions. However we simply as- 
sume that the probability distribution should be factorized as individually- 
independent one as in Eq. (5). As seen above, under this approximation, 
we can characterize the system by use of the magnetization m. In other 
words, the surrounding spins adjacent to a particular spin being replaced 
by a uniform variable m, we convert the many-body problem into a simple 
one-body system. Let us rewrite the Hamiltonian by the replacement of the 
surrounding spins Sj = m to Si as 

T N 

H MF (S i ) = -^J2 S " ( 21 ) 
1 i=i 

where z is the coordination number (usually in the case of the hyper cubic 
system, z = 2d, where d is the dimension). Similar calculations to the 
previous example give rise to the final expression of the probability for m 

as 

P/v(m) cx J dm cxp |iVmm + N log ^2 cosh (^^-m — | • (22) 
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m 

Fig. 2. Behavior of f(m,rh*) for z = 6, while we set J = 1, with /3 = 0.0 (solid), 
0.1 (dashed curve), and 0.25 (dotted curve). 



The saddle point equation is 

Oj I J ZTTI \ 

-jp- = -> m = tanh f ,3— ml. (23) 

Thus we find the probability of the magnetization for the interacting system 



Pjv(m) oc exp <N ( m 2 — mtanh 1 m — log = ) 1 

L V 2 VI -to 2 / J 



(24) 



We describe the behavior of the free energy as in Fig. 2 Readers can find 
two maximizers in the low-temperature region. Let us evaluate explicitly 
this unexpected result. The most probable state is given by the maximizcr 
of the following self-consistent equation 

m = tanh(/3 Jzm). (25) 

The explicit value of the magnetization over all temperature is shown in 
Fig. 3. Beyond the special point f3 = 1/zJ, the magnetization suddenly 
take non-zero value despite of absence of magnetic field. This is the sponta- 
neous symmetry breaking, in other terms, the phase transition. This is also 
a peculiar property in many-body systems. The mean-field analysis gives 
several meaningful properties including the phase transition on many-body 
systems as above. The results by the mean-field analysis can be validated in 
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Fig. 3. Magnetization given by the mean-field analysis for the Ising model z — 6. 
Here we set J — 1. 

infinite dimensions and reflect on the behavior in such higher dimensions. 
Therefore we must take care of their applicability to the actual behavior of 
the many-body system in finite dimensions. 

2.5. Phase transition 

Let us more qualitatively discuss the phase transition through the results 
obtained by mean-field analysis. Without the magnetic field, the spins do 
not have any tendency of the direction, that is m « 0. Indeed the most 
probable state is given by m — 0. This solution describes the paramagnetic 
phase of the magnetic material without any magnetization. However if you 
expand the self-consistent equation under the assumption that m should be 
small as 



We find the special temperature f3 — 1/zJ, at which the magnetization 
m can take a non-zero value. This is a signature of the phase transition 
from the paramagnetic state into the ferromagnetic state. We expand the 
self-consistent equation up to third order of m, and thus obtain 



m — (3Jzm. 



(26) 




(27) 
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This equation indeed gives the ferromagnetic solutions with non-zero mag- 
netization as m = ±^/3(l — f3Jz)/(f3Jz), which vanish at (3 = 1/zJ. This 
special temperature is called the critical point. In addition, two phases, 
paramagnetic and ferromagnetic phases are characterized by the behavior 
of the magnetization. Such a quantity, which can distinguish the phase, is 
called as the order parameter. Could you find the fact that you have two 
ferromagnetic solutions of the magnetization? It implies that the two max- 
imizers, say two probable states, are allowed to appear in Eq. (24). Then, 
which solution should be chosen and will appear naturally in experience? 
It depends on the dynamics rule and the initial conditions. If you design 
the dynamics to simulate thermal fluctuations At least, we can find the fact 
that, in the ferromagnetic phase, state with m = would not be expected 
to occur in the case with a large N. That means that two of the proba- 
ble states are separate and can not change one into another in the actual 
experience. In other words, the system has some robustness against ther- 
mal fluctuations. That will be an analogous property to protect the fragile 
quantum state later. 

2.6. Spin Glasses 

The above exemplified model have a uniform interactions between the Ising 
spins. The theoretical model of spin glasses is usually not the case. The 
simplest model, the Edwards- Anderson model, is defined by the following 
Hamiltonian 

H = — JjjSjSj, (28) 

where is the disordered interactions assumed to obey several type of the 
distribution functions. We often deal with two typical cases, which take a 
bimodal distribution (then called as ± J Ising model) 

P(Jij) = P 5(J, l3 - J) + (1 - p)S(J + J {j ), (29) 

and the Gaussian distribution with the average Jo and the variance J 

P(J«) = ^exp(- 5 L(J«-J ) a ). (30) 
The partition function should be dependent of the specific configuration as 

m Va}) = E n cx p (MjSiSi) • ( 31 ) 

{S.} (ij) 
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In order to evaluate the partition function and thus the free energy, we 
strive the difficult task to deal with the non-uniform interactions. Instead, 
we usually take a wise strategy to evaluate the averaged free energy based 
on the self-averaging property as, in the large-limit TV, 

^F{p-{J tJ })^^[F{p-{J l3 })], (32) 

where the square bracket denotes the average over all the combinations of 
{Jij} (configurational average) and the free energy is defined as 

-0F(fi; {Jij}) = logZGS; {Jy}). (33) 

The self-averaging property is valid for other observables, which can be 
obtained from the free energy per site. For instance, the magnetization per 
site satisfies 

I N 

i=i 

The average through the logarithmic term is still the intractable task. The 
replica method is then useful to perform the configurational average for the 
free energy. First, we evaluate the averaged power of the partition function 
to the natural number n as [Z n (j3; {Jij})} - Then we consider the analytical 
continuation of n and take the limit based on the elementary identity as 

[logZ(/3; {Jij})] = lim l^MMizl . (35) 

n->0 n 

You see the replica method everywhere in analyses on spin glasses. 

We avoid the complicated details of analytical results in spin glasses 
in order to straightforwardly understand the most important parts in this 
chapter. However we write down a few things on spin glasses related with 
our topics. The peculiar feature in spin glasses is the extraordinary slow 
relaxation toward equilibrium state in the low temperature. This is because 
the existence of many minima of the free energy on spin glasses revealed by 
the mean-field analysis. This fact implies that there are a large number of 
the most probable states. The number diverges exponentially as increase of 
the number of spins N . Unfortunately, this fact has been confirmed by the 
mean-field analysis. We have still not given the answer whether the actual 
finite-dimensional spin glasses follow the same scenario as that given by 
the mean-field analysis or not. The main reason of lack of understanding 
of finite dimensional spin glasses is absence of systematic tools to approach 
the issue. One of the exceptional methods would be the gauge theory. 
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2.7. Gauge theory 

We take the ±J Ising model as an instance of the spin glass to show the 
detailed analysis. For simplicity, we combine the strength interaction with 
the inverse temperature as K = f3J. We write the Hamiltonian of the ±J 
Ising model as 

H = -^JnjSiSj. (36) 

m 

The sign of the interactions follow the distribution function 

P(nj) = P6(nj - 1) + (1 - p)S(r ij + 1). (37) 

The partition function can be also written as Z(K; {r^ }). 

We then define the local transformation by a binary variable a% = ±1, 
called as the gauge transformation, as 7,22 

Tjj — > <7i<7jTij (38) 

Si -> ct.S,. (39) 

Notice that this transformation does not change the value of the physical 
quantity given by the average over and Si since it alters only the order of 
the summations. The Hamiltonian can not change its form after the gauge 
transformation since the right-hand side is evaluated as 

- ^2 jTijUiGjUiSiGjSj = H. (40) 

<«> 

As this case, if the physical quantity is invariant under the gauge transfor- 
mation (gauge invariant), we can evaluate its exact value even for finite- 
dimensional spin glasses. The key point of the analysis by the gauge trans- 
formation is on the form of the distribution function. Before performing the 
gauge transformation, the distribution function can take the following form 
as 

K T 

P ^ - 2^hV <41) 
where exp(— 2K P ) = (l—p)/p. The gauge transformation changes this form, 
in a different way from the Hamiltonian, as 

Let us evaluate the internal energy by aid of the gauge transformation 
here. The thermal average of the observables O is defined as, by use of the 
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canonical distribution 



Then the thermal average of the Hamiltonian can be written as 



(43) 



= £ Z(K7^ H R exp{KnjSiSj) (44) 



= -J— log Z(K;{r l3 }). (45) 



{Si} ' J j ' 

A 

dk 

We can use the self-averaging property here since this is given by the deriva- 
tive of the free energy, and thus take the configurational average as 

{Tij} (ij) 

where [■ • • ]k p denotes the configurational average with K p . Then we per- 
form the gauge transformation, which does not change the value of the 
internal energy 

W^En l^r' ^H), (47) 

Therefore we here take the summation over all the possible configurations 
of {<Ji} and divide it by 2 N (the number of configurations) as 

{»i}{T«}<«> P 

We take the summation over {(7^} in advance of that over {r^} and then 
find the partition function with K p instead of K. 

{ T ij } 

Going back to Eq. (44), we can delete both of the partition functions on 
the denominator and numerator by setting K v — K as 

[{h)k]k = 2^0^ £ £ jk cxp 

{Si} {Tij} 

= -N B tanhK. (50) 

Similarly, we can evaluate the rigorous upper bound on the specific heat as 
well as the restriction on the structure of the phase diagram. The condition 
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K p = K defines the special subspace in which we can perform the exact 
analysis even for finite-dimensional spin glasses. This subspace is called 
as the Nishimori line. 7,22 On this subspace, we can reveal several rigorous 
properties on the structure of the phase diagram even for spin glasses by 
relatively simple calculations. For instance, let us consider to evaluate the 
local magnetization (Si) by the gauge transformation. The local magneti- 
zation identifies the existence of the ferromagnetic order and its definition 
is 

/C\ \ . 1 T CXpjKTijSjSj) 

(Si)K-^Si[[ . (51) 

{s t } (ij) y ,x l]S> 

Let us consider to evaluate its configurational average 

{r«}«> P 

After the gauge transformation, we obtain 

k*>«i«, =s>n ^ar' ■<<*>«. (53) 

By summing over all the possible configurations of the gauge variables and 
dividing the resulting equality by 2^, we reach 

w ^-g, 4 ( S)' ,w ^ w " (54) 

On the other hand, let us evaluate the configurational average of the product 
of the correlation functions with different temperatures as 

[^) Kp (S i)K \ Kp = E Ill^g 1 >< (55) 

On this equality, the gauge transformation only changes the sign of the 
interactions as — > Tija^a'j in the exponential function. We find the fol- 
lowing equality by summing over {cr-} and multiplying 1/2 N 

l(Si)K} Kp =[(<Ti)K p {Si)K] Kp . (56) 

Let us discuss the structure of the phase diagram of the ± J Ising model by 
use of this equality. Setting K = K p , we obtain 

[<*)k,]k 9 = [ W **]k,- (57) 
The quantity on the right-hand side is the order parameter of the spin- 
glass phase, called as the spin-glass parameter. In the spin-glass phase, the 
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directions of the spins orient are random in space. Thus m should be zero 
and it implies 



On the other hand, if the spins are frozen, their square should be non-zero. 
Thus 



Therefore, on the Nishimori line, we prove that the spin-glass phase does 
not exist since the important fact q = m is given by Eq. (57). In addition, 
if we take the absolute value of the magnetization away from the Nishimori 
line, we find 



I [{Si)K] Kp I < I [{Si) K ] Kp I x I [(Si) Kp ] Kp | < | [{Si) Kp \ Kp I- (60) 



This inequality states that the magnetization takes the largest value on the 
Nishimori line along the vertical line of K p . It implies that a special critical 
point, multicritical point, is located at the most left-hand side of the phase 
diagram as in Fig. 4. 

In the following sections, we will find several problems on spin glasses, 
which are lying between statistical mechanics and quantum information 
processing. Although it is not enough to understand the whole properties 
on spin glasses, let us first go to these fascinating parts. After showing the 
interesting connections, we will strive the specific problems related with the 
spin glasses we must solve then. 

3. Quantum error correction: surface code 

In the first part of this chapter, we describe a technique for the quantum 
error correction. In order to protect the vulnerable quantum state from 
decoherence, we consider a method for circumvention of the effects from 
decoherence. It is very important to perform the quantum information pro- 
cessing. We here show an elaborated technique by use of the property of 
the topology. Our approach, in short terms, is based on encoding a few 
"logical" qubits in a particular state, which is not disturbed directly by 
noise of the "physical" qubits. Here we call the quantum state describing 
superposition of the binary state and 1 (in terms of informatics, bits) as 
qubits This strategy is analogous to the classical counterpart, in which we 
introduce some redundancy to restore the original state. 




(58) 



?0. 



(59) 
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Fig. 4. Nishimori line and the typical phase diagram of the ± J Ising model. The vertical 
axis is the temperature, since K can be regarded as the inverse temperature. The hori- 
zontal axis expresses the density of the ferromagnetic interactions in terms of K p . The 
ferromagnetic, paramagnetic, and spin-glass phases are denoted by "F" , "P" and "SG" , 
respectively. The solid curve represent the phase boundary. The dashed curve describes 
the Nishimori line. The multicritical point is located at the most left-hand side. 



3.1. Error model 

The quantum state of the single qubit is written as 

|*) = a\O)+0\l). (61) 

Here the above coefficients follow \a\ 2 + \/3\ 2 = 1. We can write all the 
changes on the binary Hilbert space by combination of the identity operator 
and the Pauli operators X, Y and Z. 

*-(:;)■ "-(';)■ -0-0- «»> 

By this mean, we express the error by the action of the Pauli operator. The 
action of X represents a phase error, and that of Z expresses a flip error. 
In addition, Y is a multiple error as Y = iXZ . 

Assuming the error will occur in a stochastic manner, let us define a 
noise model where the qubit gets errors as 



P -> PiP + (pxXpX + pyYpY + p z ZpZ) 



(63) 
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Although we can deal with any cases of pi^XtPy, and pz, Let us start 
from the simple case with uncorrelated between the flip and phase errors 
as pi = (1 — p) 2 and px = Pz = P, while py = p 2 , where < p < 1. In this 
case, we can independently treat two of the errors. We must construct an 
ingenious procedure to recover the damaged qubits, that is the quantum 
error correction. Once an error occurs on the single qubit system, we can not 
remove the error, since we do not know the original state. The single qubits 
is too weak to save the specific information once error occurs. Therefore we 
need to prepare many qubits against errors as the first strategy. The first 
stage of the quantum error correction is thus to construct the many-body 
quantum system to store several bits of the information. That is analogous 
with the concept of redundancy in the classical information. 

3.2. Surface code 

Let us construct an array of qubits on a torus by setting qubits on each 
edge (ij) of the square lattice embedded on a torus (genus 1). Here we 
make two ways to describe the square lattice, while being unchanged of the 
location of the qubits, as in Fig. 5. The original square lattice is convenient 
to explain how to detect and remove the flip errors and the other (dual) 
is for the phase errors. The embedded qubits (say physical qubits) on all 
the edges actually suffer from disturbance by noisy environment. We do 
not directly use the quantum state of these physical qubits to encode the 
information. Instead, we make the particular quantum state, which is stable 
against direct effects due to the errors. In order to construct such a fault- 
tolerant subspace, called as the codespace, for encoding the information, 
we define two of the check operators. One is the star operator for each site 
s as 

X s = <S>(ij)esX( lj ), (64) 
and the other is the plaquette operator for each plaquette p 

Z p = ®{ij)e P z {ij)i ( 65 ) 

where the product consists of four edges adjacent to each site and each 
plaquette as depicted in Fig. 6. All of the star operators X s and the pla- 
quette operators Z p commute and are thus simultaneously diagonalizable. 
The codespace to encode the specific information consists of their simulta- 
neous +1 eigenstates of all the check operators as X S |\I/ C ) = +1|<I' C ) for any 
s and Z p \^ c ) — +l\ 1 i> c ) for any p. The site on the original square lattice is 
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Fig. 5. Square lattice (say original) and its dual on a torus. The location of the qubit 
is the center of the edge (ij) (or denoted by the black square. 
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Fig. 6. Star and plaquettc operators. 



replaced by the plaquette on the dual one, and vice versa. The star opera- 
tors are thus found to consist of the unit loops on the "dual" square lattices. 
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Fig. 7. Contractible loops C and C* and their expressions by the check operators X s 
and Z„. 



Let us consider contractible loops on the surface of the torus denoted by C 
and C* on the original and dual lattices, respectively. The action C* can 
be expressed by the product of X s as in Fig. 7. Here notice that the dou- 
ble action of the same Pauli operator becomes unity. On the other hand, 
the plaquette operators conform the unit loop on the original square lattice. 
Thus the action represented by any contractible loops on the original lattice 
C is given by the product of Z p as in Fig. 7. The action described by C and 
C* acts trivially on the codespace, since the codespace is the eigenspace 
of the star and plaquette operators. In other words, the compatible-loop 
actions by the Pauli operator Z on the original lattice and X on the dual 
one can not alter the particular quantum state encoded on the torus. 

We here show the explicit form of the encoded quantum state below. Let 
us prepare a uniform linear-combination state of all the possible contractible 
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loops with the products of X s and Zp as 

V S Si,S 2 / 

x l + ^Z p |^2 Pl 2 P2 + ... |$„), (66) 

\ P Pl>P2 / 

where |*o) is the vacuum state of the physical qubits. Then we can easily 
confirm that X s \^ c ) = |* c ) and Z p \^ c ) = |* c ). 

On the other hand, any non-contractible loop, winding around the torus, 
of the Pauli operators of X and Z can map the codespace to itself in 
a nontrivial manner, since the possible combinations of the contractible 
loops never constitute any non-contractible loops, while such an operator 
can commute with any star and plaquette operators. If we set L x L lattice 
on a torus, we have 21? physical qubits and 2(L 2 — 1) check operators. The 
remaining degrees of freedom of 2 implies existence of two non-contractible 
loops, winding around the hole of the torus L v and winding around the body 
of the torus L t as depicted in Fig. 8. Let us express these non-contractible 
loops on the torus in terms of the products of Pauli operators as 

X v = H X m (67) 

(ij)eL v 

Xt= II ( 68 ) 

{ij)eL t 

Zv= n %?) ( 69 ) 
z t = n ( ? °) 

They are termed as logical operators. Here we use the asterisk denoting 
the non-contractible loop on the dual lattice. The combinations of non- 
contractible loops yield 2 4 = 16 different homology classes embedded in 
the original and dual square lattices on the single torus. The elementary 
manipulation confirms that the logical operators can form Pauli algebra of 
two effective qubits encoded in the topological degrees of freedom on the 
torus as [Z v , Z t ] — [X v ,X t ] — 0, while X v Z t = —Z t X v and X t Z v = —Z v X t . 
Thus, by use of these algebras by the non-contractible loops on the torus, 
we can prepare a 4 x 4 Hilbert space, say two logical qubits, explicitly given 
by |*fci,fe,,fcs,fc4> = Z^Z^X^X^ C ), where h is the number taking 
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z v 



x t 



Fig. 8. Non-contractible loops L v and Lt and the logical operators. 



and 1 to distinguish the basis of the logical qubits. The logical qubits is in 
general written as 



k\ ,&2 



(71) 



where A^kaMM is the coefficient following ^^MMM l^MMMl 2 = 1 - 
We can indeed prepare the particular quantum state from many redun- 
dant qubits as shown above. Our next interest should be how this quantum 
subspace is stable against the errors on the physical qubits. 



3.2.1. Check operators and error syndrome 

The effect coming from the decoherence, error, can occur everywhere on 
the torus. Let us assume that the error can be individually independently 
generated on each physical qubit on the torus through the channel model 
defined in Eq. (63). The errors Z/ij) and Xufi can be described as, in 
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Fig. 9. Detection of the error by the check operator. The left panel shows the original 
error chain E on the edge 1. The right panel describes the outputs of the star operators, 
which describe the endpoints of the error dE. The white circles denote nontrivial outputs 
— 1 and the black ones represent +1. 

general, non-closed error chains E and E* on both of the original and dual 
lattices. We should detect the location of the errors to circumvent their 
effects. The star and plaquette operators can also play a roll to generate 
the partial information of the errors on the physical qubits. Thus these are 
often called the check operators. The endpoints of the error chains dE and 
dE* can be detected by applications of star and plaquette operators due to 
anti-commutation of error with adjacent operators as follows. For instance, 
let us consider the local error by the action of Z on the original square 
lattice as in Fig. 9. The action of the error can be written as Z\\^). If we 
apply the star operator around Z\ on the square lattice, we obtain a pair 
of the outputs with the negative value — 1 as 

X x X^X A Z x \^l) = -ZxXxXiXiX^) = -Zi|*> (72) 
X 1 X 5 X 6 X r Z 1 \^f) = -ZxXrX^X^) = -Zi|tf>, (73) 

where we use anti commutation as Z X X X = —X\Z\. On the other hand, 
the flip error by X can be detected by the plaquette operator on the dual 
lattice. At least, we can identify the location of the endpoints on the error 
chains denoted as dE and dE*. 

We do not know the exact shape of the error chains E and E* . We must 
remove the effect of the errors on the torus only by use of the outputs of 
the check operators, error syndrome, dE and dE*. Thus we must infer the 
original state from the damaged one with negative outputs on several sites 
and plaquettes. To remove the damage signaled by the negative eigenval- 
ues, we consider to apply additional actions of the product of the Pauli 
operators. This procedure is to simply connect two endpoints of the error 
chains we can know. The additional actions with the original error chains 
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can conform closed loops. The resulting effects are trivial on the codespace. 
However, if nontrivial, we can no longer restore the original quantum state. 

In order to make the resulting effect trivial, we have to choose additional 
actions while inferring the same homology class as the logical operators 
encoding the original information only with the knowledge of the endpoints. 
The best way (optimal procedure) is to find out the most probable homology 
class. It reads, if we define P(E, E*\dE, dE*) as the probability of the 
inferred homology class E and E* conditioned on dE and dE* , 

maxE,E*P(E,E*\dE,dE*), (74) 

where E and E* denote the possible homology class of the error chains. 
Observant readers begin to recognize an existence of the limitation of the 
above procedure. First, increase of the error probability p might allow the 
possibility of the non-contractible error chains wounding the torus. If we 
prepare sufficient large torus, this is not expected to occur as <~ p L — > 0. 
As this instance, the system size L becomes larger, tolerance of the surface 
code against the error can be enhanced. However, if p is not small, we can 
not always determine the same homology class with the error chains among 
several candidates, since we tend to infer the wrong homology class due to 
the existence of long and many error chains. The connections among the 
original error chains and the additional actions would conform a long loop 
wounding the torus if p exceeds a threshold. Such an accuracy threshold 
actually exists and can be determined by special technique developed in 
the theory on the phase transition in statistical mechanics. In practice, the 
ideal error correction by inferring the most probable homology class of the 
error chains should be too harmful due to computational costs. We need to 
take some approximate technique in a moderate time. However, if we can, 
it is important to estimate the optimal error threshold in order to reveal 
the theoretical limitation of the surface code. 

In the following, we construct the statistical-mechanical model to pave 
the way to analytically estimate the optimal accuracy threshold. 

3.2.2. Probability of error chains 

Notice that our task is to elucidate the macroscopic property in the error 
correcting procedure, not actual many-body quantum system. The outputs 
of the measurements of the check operators and the procedure of the ap- 
plication of additional actions based on the inferred homology class are 
classical. Although we deal with the problem on the quantum error correc- 
tion, all the methods we employ here may not be quantum. 
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E + E' = C 



Fig. 10. Reasonable connections E' inferred from dE. For transparency, the square 
lattice is depicted by the dotted lines. 



The error on each qubits occur in the probabilistic manner. Therefore 
let us construct a probabilistic model of the quantum error correction on 
the torus. In order to deal with the asymptotic behavior of the error pattern 
and inferred chains in the infinite- number limit, we use the most suited tech- 
nique, statistical mechanics. Statistical mechanics starts from constructing 
a probabilistic model in the microscopic level. Following the concept of the 
large deviation principle, we take an infinite number limit of the system 
size in order to predict the deterministic macroscopic behavior. 

Now we take the simple case with uncorrelated errors. We can separately 
deal with the flip and phase errors. For simplicity, hereafter we consider only 
the phase errors given by the actions of Z on the original square lattice. 

From the knowledge of endpoints dE, error syndrome, we have to infer 
the most likely homology class of the error chains, while considering any 
reasonable choices. 

For instance, let us consider the reasonable connections E' between each 
pair of the endpoints dE as in Fig. 10. As a result, we can make contractible 
loops denoted by C on the square lattice. Then the reasonable chains E' 
inferred from dE are in an equivalent class with the error chains. However 
there are several cases to infer non-equivalent class with the error chains as 
in Fig. 11. 

By use of the Bayes theorem, we consider to evaluate the probability 
to infer the equivalent homology class with the original error chains condi- 
tioned on those endpoints as follows 

P(dE\E)P(E) 



P{E\dE) = 



(75) 



where Dk 1 ,k 2 ,k 3 ,k l denotes the logical operator. The equivalent class is 
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E E'.dE E + E 1 = C + A,,i,o,o 

Fig. 11. Failure of inference of the equivalent class with E. The wounding chain by the 
product of Z yields Zt denoted by (ki, ki, k%, ki) = (0, 1, 0, 0). 




given simply by -Do, 0,0,0 = 4>- 11 we infer the different class denoted by 
(k 1 ,k2,k 3 ,k 4 ) ^ (0,0,0,0), the error correction fails. We assume that the 
prior probability for the homology classes is uniform, namely P(dE) = 
1/16, since we encode various combinations of the different homology 
classes. Our task is to consider P(dE\E), which is the probability to gen- 
erate the endpoints of the specific error chains. Let us first evaluate the 
probability for the specific error chains. The error chains are generated 
following the independently identical distribution as 



m v 7 



where tE represents the error chains and takes ±1 {tE < 0, when (ij) € E). 
Here we introduce the expression of the exponential form as in the case for 
the spin glass and obtain 

n exp(^) 



(ij) 



I' 



Notice that K p is minus of the original definition in the spin glass. Then, 
in order to evaluate P(dE\E), we sum over all the possible configurations 
E' = E + C of the error chains sharing the endpoints. We reach 

P{dE\E) « ^U^Kpr^rf/), (78) 

c (ij) 
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where the summation is taken over all the possibilities of C and the product 
is over all the edges. We use the same indicators as r^j for C. 

This probability tells us how to proceed the error correction. We simply 
perform the connection of the endpoints in the stochastic manner. If we 
know p in advance, it is better to set the same value. Otherwise we have 

P K (dE\E) cx Y.W^MKr^). (79) 

c (ij) 

The parameter K stands for the importance/preference to choose the rea- 
sonable choice in the connection of the endpoints. For instance, in the case 
with K — > oo, E' with the minimum length are preferred. The performance 
of the error correction depends on the value of K. 

The loop constraints II (ij) T ij = 1 allow to use another expression by 
the Ising variables tF.n = SiSj. By use of this expression, we can find that 
P{E) is written by the partition function of the ± J Ising model as 

Pk(8E\E) cx Z(K; { Tij }) = £J] > ( 80 ) 

{s>} (ij) 

where is the sign of the random coupling in context of spin glasses. 
Each of the random couplings follows the distribution function of the er- 
ror chains P{E). Notice that, in this case, p denotes the density of the 
antifcrromagnetic interactions. Therefore the problem to identify the error 
threshold of the surface code is reduced into evaluating the quantity re- 
lated with the partition function of the simple spin-glass model, the ± J 
Ising model on the square lattice. If we increase p, the spin-glass system 
will be not rigid against thermal fluctuations, since the ferromagnetic order 
decays. Similarly, in the surface code, increase of p is expected to lead the 
error correction to be unfeasible. These analogous implies the existence of a 
fascinating relationship between the phase transition in the spin glass and 
the error threshold in the surface code. 

In order to more clearly see this connection between statistical mechan- 
ics and quantum error correction, we show a schematic picture of the phase 
diagram for the ±J Ising model on the square lattice as in Fig. 12. In the 
low-temperature region with a relatively small p, the ferromagnetic ordered 
phase can be observed. It means that the order of the Ising spins exists and 
suppresses the fluctuation of the domain walls (boundaries between different 
signed spins). Let us consider to evaluate Pk{E\8E) in the ferromagnetic 
phase. We rewrite the probability as, in terms of the free energy of the spin 
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Fig. 12. Phase diagram of the ±J Ising model. The figure is depicted by the same 
symbols in Fig. 4. In the two-dimensional case, the absence of the spin-glass phase has 
been partially proved and supported by numerical simulations. 23 



glass, 

= 1 + E| „e*p(W W .) - (81) 
Here the difference of the free energy AF^^MM is given as 

&F klMMM = F(K; { Tij }) - F(K; {r t3 T^ k ^}), (82) 
where we use K instead of j3 as 

F(K;{ Tl] r^- k ^}) = log '/. : K: {","/,'' : ; ! )• (83) 

In context of statistical mechanics, the logical operator Dh^MMM i s i n ~ 
terpreted as the induction of the antiferromagnetic interactions along the 
non-contractible loop as {TijT i - kl ' k2 ' k3 ' ki } . In the ferromagnetic phase, such 
a long boundary of antiferromagnetic interactions yields a crevice with the 
same scale as the inducted loop in the ordered spins. Therefore the free 
energy difference should become ~ O(L). As a result, the probability to 
infer the equivalent homology class with the original error chains is asymp- 
totically, in the limit of L — > oo 

P K (E\dE) ->■ 1. (84) 



March 2, 2013 16:21 



WSPC - Proceedings Trim Size: 9in x 6in O2'0hzeki 



31 



This means that the error correction is feasible. 

On the other hand, in a high temperature and for not a small p, namely 
the paramagnetic phase, the order of the spins is destroyed. The vast num- 
ber of islands of the different-signed spins can exist and the boundaries 
fluctuate. The free-energy difference becomes zero, and thus we reach. 



This means that the error correction is infeasible. Therefore the locations of 
the critical points of the spin-glass model identify the error thresholds of the 
surface code. In particular, the optimal threshold is located on the special 
critical point along the Nishimori line, namely multicritical point, 7,22 since 
it is at the most right side in the phase diagram. In other words, if we know 
the precise value of p, it is easier to correct the errors in the logical qubits. 

3.3. Analyses on accuracy thresholds for surface code 

In general, the problem on finite-dimensional spin glasses is intractable. 
However recent development in theory of this realm enables us to estimate 
a precise location of the critical point in several spin glasses. In the following 
section, we introduce such a specialized theory in detail. 

3.3.1. Duality analysis: simple case 

The situation that no systematic analytical methods attacking the prob- 
lems on the critical phenomena in finite-dimensional spin glasses have been 
changed since a recent development in the spin-glass theory. It enables us 
to estimate the precise value of the special critical point especially on the 
Nishimori line, which corresponds to the optimal error threshold. 24-27 The 
method as shown below is based on the duality, which can identify the lo- 
cation of the critical point especially on two-dimensional spin systems. 28 ' 29 
Let us review the simple case of the Ising model. 

The duality is a symmetry argument by considering the low and high- 
temperature expansions of the partition function. As a result, we can obtain 
a simple relation between two different temperatures through the partition 
function as 



P K {E\dE) -> 1/16. 



(85) 



Z{K) = AZ(K*), 



(86) 



where K* denotes the dual coupling constant related with the original one 
K . The coefficient A will be an important quantity below. 
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We here deal with the case for the Ising model on the square lattice, 
whose partition function is given as, through the Hamiltonian (20), 

Z(if) = ^H ex P(^ 5 i). ( 87 ) 

where the product is taken over all the nearest neighboring pairs on the 
square lattice. The original formulation of the duality analysis is based 
on a relatively painful calculation. 28 We here show a much simpler version 
given by a simple Fourier transformation for the local part of the Boltzmann 
factor. 29 

We define the edge Boltzmann factor as 

x^K) = cxp(K cos7r</>), (88) 

where we use another binary variable as <fr — and 1, instead of the Ising 
spin. In addition, we apply the binary Fourier transformation to this quan- 
tity, called as the dual edge Boltzmann factor, 

x* k {K) = Y] exp(i7rfcr/>) = (c K + c~ K cos?rfc) . (89) 
^ 2 0=0,1 ^ 2 

Then we find that the partition function can be written in two ways by 
both of the edge Boltzmann factors. First, the partition function is simply 
expressed by the original Boltzmann factor as 

w=eii^W' ( 9 °) 

where <fnj — <fn — 4>j. The difference between (j>i and <j)j is taken in order 
from left to right and from top to bottom. Inserting the inverse Fourier 
transformation of the dual edge Boltzmann factor yields another expression 
of the partition function 

Z W = [-7=) EE H^iKy^-^. (9i) 

{4>i} {ha} (ij) 



Here let us perform the summation over {&} by considering the adjacent 
edges. We concentrate on the particular site and then find that the expo- 
nential term can be factorized and the summation can be taken indepen- 
dently as 

c i(-fcio-fc2o+fco3+fc O4 )0o = 2<5(fc 10 + k 20 - k 03 - k 04 = (mod 2)). (92) 

00=0,1 

We describe adjacent edges to the particular site as in Fig. 13. In order 
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k 



20 



k 



10 



0; 



^04 
&03 



Fig. 13. Description of the terms appearing in Eq. (92) 

to remove the terms of the Kronecker's delta functions, we employ another 
binary variable as 



k20 = 

ko3 = 
k 0i = 



(93) 
(94) 
(95) 
(96) 



We set these new variables {</>*} on each plaquette on the square lattice, 
namely each site on the dual one. The resultant expression of the partition 
function is 



Wife) 

This fact leads us to the double expressions with the 



where 4>* 3 = <p* - 4>* 
different arguments as 

Z(x (K), Xl (K)) = Z{x* (K),xl(K)). (98) 

To reduce the number of the arguments, we normalize the partition function 
by the principal edge Boltzmann factors xo(K) and Xq(K). 

{x (K)} Nb z{u x {K)) = {xq(K)} Nb z(ul(K)), (99) 

where z is the normalized partition function z(u\) = Z/{xo(K)} Nn and 
z(u\) = Z/{x* (K)} Nb . We explicitly obtain ui(K) = xi(K)/x (K) = 
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exp(-2K) and u*{K) = x*(K)/xq(K) = t arm if. It reads, if we set the 
dual coupling as exp(— 2K*~) = tanh(AT), 

{x (K)} Nb z{ Ul {K)) = {x*(K)} Nb z( Ul (K*)). (100) 

Therefore we find that the partition function for the Ising model has symme- 
try between two different temperatures. This is the duality. The coefficient 
A is also obtained as 

A =(5M)"* = {7i (1+exp( - 2jr)) }'" > ' (101) 

The well known duality relation exp(— 2K*) = tanhif can identify the 
exact location of the critical point by equating K = K* . The critical point 
is estimated as K c = ln(l + V2)/2 = 0.440686794 ••• . Interestingly, the 
coefficient A becomes unity at the critical point. We use this property as a 
priori assumption for the analysis on the critical point for spin glasses. 



3.3.2. Duality analysis: spin glass 

The replica method, which is often used in theoretical studies on spin 
glasses, allows to generalize the duality analysis to spin glasses. 24,25 Let us 
consider the duality for the replicated partition function as [Z n (K; {nj}}. 
In this case, the multiple binary Fourier transformation defines the dual 
edge Boltzmann factor as 

4 u k 2 ,-,k n (K p ,K)=(^-) ]T x^ 2r .., (t>n (K p ,K)e i ^ k ^, 

(102) 

where the original edge Boltzmann factor is given by the configurational 
average of the n-replicated ± J Ising model, 

x^,<h,-,*n{K p ,K) = {e- K " +K cos "** + e K "~ K ^ cos **« } . 

(103) 

Notice that the definition of p is the density of the antiferromagnetic inter- 
actions here to analyze the accuracy threshold for the surface code. They 
leads us to the double expression of the replicated partition function as, in 
a similar way to the above simple case, 

{x (K p ,K)} N °z( Ul (K p ,K),u 2 (K p ,K),---) 

= {x*(K p , K)} N -z(ul(K p , K), u* 2 (K p , K),--- ), (104) 

where the subscript of uu and u^, stands for the number of anti-parallel pair 
among n replicas on each edge. We restrict ourselves to the case to analyze 
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the location of the multicritical point, namely K — K p on the Nishimori 
line. 

Unfortunately we cannot replace ul(K,K) by Uk{K* ,K*) even after 
normalization as the above case for the Ising model, since the replicated 
partition function is multivariable. Nevertheless we can estimate the loca- 
tion of the multicritical point even without the ordinary procedure of the 
duality. We here put a priori assumption that xq(K, K) = Xq(K, K) at the 
critical point, implying the coefficient of the double expressions of the repli- 
cated partition function should be unity. According to this assumption, we 
take the limit n — > of the equation along the replica method and thus 
obtain 

-plogp- (1 -p)log(l-p) = 7; log 2. (105) 

The solution is p c = 0.1100 • • • . We thus conclude that the accuracy thresh- 
old is estimated as p c = 0.1100- • • . Strictly speaking, this analysis is not 
exact. In practical, the precision of the above result has been confirmed to 
be satisfiable by supports from numerical simulations. In addition to such 
numerical validations, the following theoretical refinement of the method 
to identify the multicritical point have been considered. 

3.3.3. Duality analysis with real-space renormalization 

As above mentioned, the ordinary duality analysis hampers since the repli- 
cated partition function was multivariable. Thus we rely on the assumption 
that Xo(K, K) = Xq(K, K) would give the location of the critical point. We 
here sketch the relationship between the double expressions of the repli- 
cated partition function as the curves of the relative Boltzmann factors 
Uk{K,K) and u* k {K,K) on the two-dimensional space for simplicity, al- 
though those are correctly in a hyper space as in Fig. 14. The thick curve 
denotes the relative Boltzmann factor Uk (K, K) and the dashed one rep- 
resents the dual one ul(K 7 K). When the temperature increases, the rep- 
resentative point of the replicated partition function moves from P (the 
high-temperature limit) to F (the low temperature limit) along the thick 
curve. On the other hand, the dual representative point goes inversely along 
the dashed curve. These features have been shown rigorously and imply the 
existence of the duality relation for the temperature 30 . If two curves become 
completely coincident with each other, we then obtain a relation implying 
u%{K,K) = iik(K*, K*). Solving this relation, we can obtain the duality 
relation for different temperatures as the case of the Ising model. In spin 
glasses, they do not overlap. 
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P Ui 

Fig. 14. Relative Boltzmann factors of the replicated partition function (Projection onto 
two dimensions). The high-temperature limit is given by P, while the low-temperature 
limit is F. The solid curve represents the original relative Boltzmann factor. The dashed 
curve denotes the dual relative Boltzmann factor. 

We import another piece of the theories in statistical mechanics, the 
real-space renormalization group analysis. Most of the problems on sta- 
tistical mechanics are tractable since the degrees of freedom are highly 
correlated. Often we employ a trick to map the original problem into much 
simpler problem with a recursive structure by use of some approximation. 
As in Fig. 15, we trace over a part of spins on the square lattice. Repeating 
this procedure while omitting the generated multi-body interactions (ap- 
proximation), we construct the renormalization group in the form of the 
recursion 

K n = R(K n _i), (106) 

where K n is the renormalized coupling constant after n steps. By use of 
the renormalization group analysis, we describe flow of the renormalized 
coupling constant in K space. The flow usually terminates two fixed points 
representing the ordered and disordered phase, while being divided by an 
unstable fixed point, which is the critical point. The precision to represent 
the original behavior in the model is dependent on the approximation in 
construction of the renormalization group. We generalize this procedure 
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Fig. 15. Real-space renormalization. In this instance, a limited number of spins are 
summed, and the remaining ones construct another square lattice with multiple-body 
interactions. 




Fig. 16. Relative Boltzmann factors of the replicated partition function and their renor- 
malization flow (Projection onto two dimensions). The dotted arrows depict the renor- 
malization flow starting from several points on both of the relative Boltzmann factors. 
The internal curves denote the renormalized relative Boltzmann factors. The bold curve 
describes the sufficiently renormalized relative Boltzmann factors with the unstable fixed 
point as the case with a single variable in the partition function. 

into the relative Boltzmann factors {uk(K, K )} and {u%(K, K)}. The flow 
can be depicted in the hyper space of {zik(K, K)} as in Fig. 16. Similarly 
to the above simple case, the renormalization flow goes toward two fixed 
points P and F. The unstable fixed point C would be located between the 
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Fig. 17. Clusters for duality with renormalization. The white circles denote edge spins 
whose configuration is represented by k, and the black ones are to be summed. 

original and dual relative Boltzmann factors, since they do not overlap but 
have a common critical point. The renormalization flow starts from both 
of the relative Boltzmann factors and once moves to the unfixed point. Af- 
ter several renormalization steps, the flow goes to two fixed points. That 
means that, if we conduct one step of the renormalization for each point 
on the original and dual relative Boltzmann factors, we observe both of 
the renormalized curves get close to each other as in Fig. 16. The sufficient 
steps of the renormalization makes both of the curves into a common thick 
curve capturing the unfixed point, namely the multicritical point. As far as 
possible, we desire to estimate the location of the critical point precisely. 
Then, let us consider to trace the partial spins without any approximations. 
For example, in the case on the square lattice, we define the cluster Boltz- 
mann factor X]? (K, K) , where the subscript k denotes the configuration 
of the edge (white-colored) spins and s expresses the size of the cluster in 
Fig. 17. When k = 0, all the edge spins are in up directions. We sum over 
the internal (black-colored) spins in order to perform the renormalization 
without approximation. It is difficult to perform the sufficient step of the 
renormalization without any approximation. Remember that we can esti- 
mate the relatively correct value of the multicritical point even without the 
renormalization. Therefore we propose a systematic way to improve pre- 
cision of the location of the multicritical point. We employ the following 
equation to estimate the location of the multicritical point, 

x ( Q s) (K,K) = X ; {s) (K,K). (107) 
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The equality for s — (edge) reproduces the case without renormalization 
as pc°^ = 0.1100. 24,25 If we increase the size of the used cluster, we can sys- 
tematically approach the exact solution for the location of the multicritical 
point of the ± J Ising model as p ( c ] = 0.1093 and p¥ ] = 0.1092. 

If we remove the condition of the Nishimori line K p — K, we can de- 
scribe the phase boundary for spin glasses by 

x { s) (K p ,K) = x* {s) (K p ,K). (108) 

By use of this equation, we can obtain the precise results for the phase 
boundary in the higher temperature region than the Nishimori line and the 
whole phase boundary of the diluted Ising model. 26 ' 27 

3.3.4. Other cases 

Not only the surface code, several quantum error correcting codes are found 
to possess the connection with the spin-glass models. Although we omit 
their detailed explanations, we look over the recent results in short below. 
We restrict ourselves to the case in which the duality analysis have predicted 
the accuracy threshold. 

First, we simply mention the cases by other arrangement of physical 
qubits in the surface code than the square lattice, say triangular and hexag- 
onal lattices. The accuracy threshold for both of the lattices can be given 
as p c — 0.164 • • • (triangular) and p c — 0.067 • • • (hexagonal). 27 Recently, 
another type of the surface code with more computational capability have 
been developed, color code. Also in the case of the color code, we prepare 
the arrangement of the physical qubits on each unit triangle on the trian- 
gular lattice or Union- Jack lattice as in Fig. 18. 31 ' 32 The color code on the 
Union- Jack lattice implements the whole Clifford group of unitary gates 
generated by the Hadamard gate, the {7r/8} gate, and the controlled-NOT 
gate, although that on the triangular lattice can not employ the {tt/8} gate. 
It means that both of the color codes have the wider computational capabil- 
ity than that of the Pauli group. The corresponding statistical mechanical 
model has three-body interactions differently from the above surface code. 
The quenched random interaction then represents the error arising on the 
physical qubits on the unit triangles. Here the duality analysis can esti- 
mate the accuracy thresholds for both of the color codes as p c = 0.1096 — 8 
(triangular) and p c = 0.1092 — 3. 33 Interestingly, the advantage of the com- 
putational capability does not spoil the robustness of the error correction. 

In addition, the practical errors in implementation of the quantum sys- 
tem are not only on the computational basis. Lack of the physical qubits is 
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Fig. 18. Color code on the triangular and Union-Jack lattices. 

also likely to occur. Recently, Stace et al. have proposed a modified scheme 
of the surface code to protect the logical quantum state from both type of 
the errors, namely from decoherence and lack of the physical qubits. 34 ' 35 
Then the corresponding statistical mechanical model becomes the diluted 
version of the ± J Ising model. We can identify the locations of the optimal 
thresholds depending on the ratio of the loss of the physical qubits q as 
in Fig. 19. 36 The detailed values are shown in Table 1. As shown in Table 
1 and Fig. 19, the decay of the robustness, depending on increase of q, of 
the surface code against the flip/phase errors can be observed. This behav- 
ior can be interpreted as the decay of the ferromagnetic order due to the 
dilution of the interactions. 



Table 1. Optimal thresholds given by the duality 
analyses with s = Q, 1, and 2 clusters for the uncor- 
rected case. 



1 


Pc (s = 0) 


Pc (s = 1) 


Pc (s = 2) 


0.00 


0.11003 


0.10928 


0.10918 


0.10 


0.09240 


0.09196 


0.09189 


0.20 


0.07245 


0.07235 


0.07233 


0.30 


0.04984 


0.05004 


0.05009 


0.40 


0.02462 


0.02492 


0.02500 


0.45 


0.01155 


0.01174 


0.01179 



3.3.5. Depolarizing channel 

Before closing this section, we add another example on the quantum error 
correction related with statistical mechanics. We discussed the case without 
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Fig. 19. Decay of the robustness due to loss of the physical qubits. The solid curve 
expresses the results by the duality with s = and the dashed curve denotes those with 
s = 1, while both of them are mostly overlapping in this scale. 



correlation of both of the flip and phase errors for simplicity. It is straight- 
forward to generalize the above procedure for the depolarizing channel case 
such that the error probabilities are equalized as px — Py = Pz = p/3- In 
this case, the corresponding statistical mechanical model becomes a little 
bit different one from the ± J Ising model. We can no longer independently 
deal with both of the errors. Let us first evaluate the probability for the spe- 
cific error chains. The error chains are generated following the distribution 
function as 



j *J I *J 

4 ~r 4 



(a) 



p 



3(1 -p) 



T E T E* T E T E* 

LL J — 12 I — Li 12 

4 4 4 



(109) 



We again use the expression of the exponential form as in the case for the 
spin glass and obtain 

P(E, E*) cx ]J exp(i<r d r| + K d rf* + K^t?*). (110) 



Notice that exp(— 4Kd) = 3(1 — p)/p. Similarly, the problem is to iden- 
tify the singularity of the following the probability through the partition 



March 2, 2013 16:21 



WSPC - Proceedings Trim Size: 9in x 6in O2'0hzeki 



42 

function of a spin glass 

P K (dE,dE*\E,E*) 

K E II ex P(^^^ + K ^*S*S* + Kr^Tf/SiSjStS*), 

Si,S? (ij) 

(111) 

where the quenched random coupling obeys the distribution function 
P(E,E*) . This model is regarded as a spin-glass version of the 4-state 
Potts model. 

The duality analysis can be performed for this statistical-mechanical 
model in order to estimate the optimal thresholds. 37 In addition, we can 
apply the modified scheme to circumvent the effects due to loss of the phys- 
ical qubits proposed by Stace et al. 36 We list the results for the depolarizing 
channel in Table 2. 



Tabic 2. Results for the corre- 
lated case (depolarizing channel) 
Px = Pz = PY = p/3. 



<3 


Pc (s = 0) 


Pc (s = 1) 


0.00 


0.18929 


0.18886 


0.10 


0.16025 


0.15985 


0.20 


0.12690 


0.12656 


0.30 


0.08844 


0.08819 


0.40 


0.04454 


0.04440 


0.45 


0.02121 


0.02114 



4. Quantum annealing and beyond 

The second part of this chapter, we introduce another application of sta- 
tistical mechanics in context of quantum computation. As in the other 
chapter, by use of the analogy of simulated annealing (SA) developed in 
statistical mechanics, a generic algorithm intended for solving optimization 
problems by use of quantum nature, quantum annealing (QA), has been 
studies since its proposal. In SA, we make use of thermal, classical, fluc- 
tuations to employ a stochastic search for the desired lowest-energy state, 
which corresponds to the optimal solution, by allowing the system to hop 
from state to state over intermediate energy barriers. In QA, by contrast, 
we introduce non-commutative operators as artificial degrees of freedom of 
quantum nature in order to induce quantum fluctuations. The most typical 
procedure of QA is performed by adiabatic control of quantum fluctuations, 



March 2, 2013 16:21 



WSPC - Proceedings Trim Size: 9in x 6in O2'0hzeki 



43 

say quantum adiabatic computation (QAC). The adiabatic theorem guar- 
antees that an sufficient slow time evolution would let the system closely 
follow the instantaneous ground state. Therefore, starting from an initial 
trivial ground state, we can reach a nontrivial one by slowly controlling 
quantum fluctuations. However, in a practical sense, we desired to perform 
the procedure to hnd the answer of the optimization problem as fast as pos- 
sible. The adiabatic control to efficiently obtain the optimal solution takes 
a characteristic time related with the energy gap between the first excited 
state and the lowest one during its procedure. The problems with closure 
of the energy gap are difficult to be solved in a moderate time by QAC. 
Unfortunately, as far as we know, the typical difficult problems in classical 
computation involves the exponential closure of the energy gap by increase 
of the problem size. 

In this chapter, we show several trials to overcome the bottleneck of 
quantum adiabatic computation by use of statistical mechanics. We use 
formal similarities between SA and QAC in order to import several devel- 
opments in statistical mechanics. In both methods, we have to drive the 
system slowly and carefully to control the strengths of thermal or quan- 
tum fluctuations. The idea behind quantum adiabatic computation is to 
keep the system close to the instantaneous ground state of a quantum sys- 
tem. This is analogous to the protocol of SA, in which one tries to make 
the system keep quasi-equilibrium state. It is indeed possible to make this 
analogy more transparent by a precise formulation, from which a fascinat- 
ing technique to be expected to improve the performance of QAC can be 
proposed. 

4.1. Quantum adiabatic computation: short review 

Let us briefly review the procedure of the most typical QA, namely QAC, 
before importing theory of statistical mechanics to quantum computation. 
In QA, we introduce a non-commutative operator to drive the system by 
quantum nature as 

H(t) = f(t)H + (l-f(t))H 1 , (112) 

where Hq is the classical Hamiltonian consisting of diagonal elements, 
which express the cost function. Here f(t) is assumed to be a monoton- 
ically increasing function satisfying /(0) = and f(T) = 1. For instance, 
f(t) = t/T, where T denotes the computation time for QAC. The quan- 
tum annealing starts from a trivial ground state of H\, which is often 
chosen to be a uniform linear combination of the computational basis as 
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1^(0)) = \a)/yN. Let us deal with only the discrete combinatorial opti- 
mization problems, which are simply termed as the optimization problem 
below. Most of the optimization problem can be expressed by the spin-glass 
Hamiltonian. In order to explain the procedure of QAC in detail, we con- 
sider to find the ground state of the simple spin-glass model as used before, 
namely the random-bond Ising model, 

H Q = -Y,JijVi°h ( 113 ) 

where the summation is taken over all the nearest-neighboring pairs of the 
Ising spins. We take the computational basis of the eigcnstatcs of the Ising 
variables to represent the instantaneous state as = | , erf, • • • , <J%)- 

The transverse-field operator is often used as quantum fluctuations for im- 
plementing QAC for the spin-glass model 

N 

ff 1 = -r 5>f, (in) 

i=l 

where To is the strength of the transverse field. The whole Hamiltonian of 
QAC (although widely used for QA) thus becomes 

N 

H(t) = /(*) £ Jirfv* + (1 - /(*)) To °i- ( 115 ) 
The ground state of the transverse-field operator Hi is trivially given by a 

i—N 

uniform linear combination as X){cr} I ")/ v 2 . For a sufficiently large T, the 
adiabatic theorem guarantees that the instantaneous state at time t, \&(t)), 
is very close to the instantaneous ground state implying |0(i)), (0(t)\^/(t)) ~ 
1, when the instantaneous ground state |0(i)) is non-degenerate. The con- 
dition for \0(t)), (0(i)|*(i))l - e 2 (e < 1) to hold is given by 



max 



<l(*)|«|0(i)) 



iA 2 (t) 



(116) 



where |l(i)) is the instantaneous first excited state, and A(t) is the energy 
gap between the ground state and first excited one. The maximum and 
minimum are evaluated between and T. In our case, since dH(t)/dt oc 
1/T, the adiabatic condition is written as 

T -^nW (117) 
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Therefore, if we desire to solve the problems involved with the exponential 
closure of the energy gap while increase of N, QAC must take extremely 
long time to find the ground state with high probability. 38,39 

4.2. Novel type of quantum annealing 

In order to overcome this problematic bottleneck of QAC, we change our 
strategy from the adiabatic control of quantum fluctuations. We demand 
an important key for noncquilibrium statistical mechanics by using a fas- 
cinating bridge between quantum computation as QAC and statistical me- 
chanics. To make the connection more clear, we show a useful technique to 
relate both of the fields. 

4.2.1. Classical quantum mapping 

Let us compare two of the procedures of SA and QAC. Both of the protocols 
are given by slow sweep of thermal and quantum fluctuations to keep the 
system trace the instantaneous stationary state, equilibrium for SA and 
ground state for QAC, respectively. In numerical implementation of SA, we 
demonstrate a stochastic dynamics driven by thermal fluctuation with the 
master equation. 



where P(a;t) is the probability with a spin configuration of {of} simply 
denoted as a at time t. Notice that a is not a state of a single spin but is a 
collection of spin states. M(a'\a; t) expresses the transition matrix following 
the conservation of probability M(a\a'; t) = 1 and the detailed balance 
condition 



where we denote the instantaneous equilibrium distribution as P cq (<r;t) = 
cxp(— f3(t)E(a; t))/Z(t) and the instantaneous energy E(a;t) is the value 
of the classical Hamiltonian H (t). Since SA consists of a dynamic control 
of the temperature, the time variable t has been written explicitly in the 
arguments of the inverse temperature and the partition function. In order 
to satisfy this condition, we often use the transition matrix with Metropolis 
rule as 




(118) 



M(a\a';t)P cq (a';t) = M(a'\a;t)P cq (a;t), 



(119) 



M(a\a';t) 



min(l, exp(-/3A£(cr|cr'; t))) 



(120) 



March 2, 2013 16:21 



WSPC - Proceedings Trim Size: 9in x 6in O2'0hzeki 



46 



where 



AE(a\a';t) = E(a;t) - E(a';t) 



(121) 



or heat-bath rule as 



M(<j\a';t) = 8x{a,a') 



exp 



(122) 



2 cosh 



(§AE(a\a>;t)) 



where 



N 



5M°') = *(2, £(1 



- 0-i<7-)). 



(123) 



On the other hand, the dynamics of QAC is governed by the Shrodinger 
equation. To look at the connection between SA and QAC, we employ 
a mapping technique of dynamics of relaxation toward equilibrium as the 
master equation into the Shrodinger equation. If we use the following special 
quantum Hamiltonian, we find it possible to simulate the dynamics of SA 
in quantum manner, 40 



H q (a'\a; t) = 5(a' , a) - e^ H °^/ 2 M(a'\a; t )c~^ H ^/ 2 . (124) 



Here we consider to gradually increase the inverse temperature following 
the spirit of SA as (3(t). This special Hamiltonian has the following state 
as its ground state, 



It is easy to confirm that the quantum expectation value of a physical 
quantity A{a) in the ground state as ("i 1 cq (t)\A(a)\^ cq (t)) coincides with 
the thermal expectation of the same quantity with (3(t). The ground state 
energy simply takes zero. This fact is shown by use of the conservation of 
the probability and the detailed-balance condition as 



On the other hand, the excited states have positive-definite eigenvalues, 
which can be confirmed by the application of the Perron-Frobenius theo- 
rem. By using the above special quantum system, we can treat a quasi- 
cquilibrium stochastic process in SA as an adiabatic dynamics as in QAC. 




(125) 




(126) 



a 
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The above formulation is a generic way of the classical-quantum map- 
ping. We demonstrate the above mapping of SA into QAC by more explicit 
instance. Let us consider an optimization problem that can be expressed as 
a classical Hamiltonian with local interaction 

H = -^Hj, (127) 

3 

where Hj involves rjj and a finite number of a^(k n eqj). Taking a familiar 
instance is a spin-glass system with nearest-neighbor interactions, 

Hj = -^2 J jk<y Z j<Jl, (128) 

kedj 

where dj denote sites adjacent to j. The following Hamiltonian is the ex- 
plicit form, which facilitates our analysis, 

fffW = -xW£(^-e^), (129) 

3 

where x(i) = e^^ p with p = maxj|ifi|. Notice that p is proportional to 
the interaction energy and is the order of O(N ) due to finiteness of the 
interaction range. Let us consider that the protocol of SA in the above 
quantum system. In the first stage with very high temperature /3(0) — >• 0, 
the quantum system (129) reduces 

^ SG (*) = -E(^ -!)■ (130) 

3 

Its ground state is the uniform linear combination of all possible states 
in the basis to diagonalize {erf}. It means that all states appear with an 
equal probability as the equilibrium distribution in high-temperature limit. 
Therefore the quantum system (129) can correctly demonstrate the initial 
condition of SA. In addition, in the limit (3(t — > T) — » oo, (129) becomes 

3 

The ground state is with the lowest value of the classical system (127), be- 
cause each Hj takes its lowest value. These observations confirm that the 
quantum system (129) indeed demonstrate the quasi-stationary dynamics in 
SA. Notice that, Interestingly, the adiabatic condition of the above special 
quantum system we used in the classical-quantum mapping can reproduce 
the condition of convergence of SA. By use of the fascinating connection be- 
tween SA and QAC, we can import several theories of statistical mechanics 
into the quantum dynamics. 
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As shown later, the collaboration of statistical mechanics with quantum 
dynamics can produce a new algorithm to solve optimization problems in 
different manners from QAC. 

4.2.2. Jarzynski equality 

Among several recent developments in statistical mechanics, we take the 
Jarzynski equality (JE) as an attempt to improve the performance of 
QAC. 41 This equality relates quantities at two different thermal equilib- 
rium states with those of nonequilibrium processes connecting these two 
states. It can also be termed as a generalization of the well-known inequal- 
ity, the second law of thermodynamics AF < (W)o->t- Here the brackets 
(• • ■ )o->t arc for the average taken over nonequilibrium processes between 
the initial (at t = 0) and final states (at t — T), which are specified only 
macroscopically and thus there can be a number of microscopic realizations. 
The Jarzynski equality is written as 42 ' 43 

(e-^)o^ = f ■ (132) 

Here the partition functions for the initial and final Hamiltonians are ex- 
pressed as Z and Zt , respectively. One of the important features is that JE 
holds independently of the pre-determined schedule of the nonequilibrium 
process. Another celebrated benefit is that JE reproduces the second law of 
thermodynamics by using the Jensen inequality. Notice that we have to take 
all fluctuations into account in evaluation of the expectation value in the 
right-hand side of JE in order to calculate the free energy difference. The 
Jarzynski equality holds formally in the case with change of temperature, 

<Oo->T = f. (133) 

when we employ the pseudo work instead of the ordinary performed work 
due to the energy difference as 

Y(a;t k ) = (0 k+1 -0 k )E(v;ti). (134) 

Here we employ discrete time expressions as to — and t n — T for simplic- 
ity. We show the simple proof of JE for the dynamics in SA. Let us consider 
a thermal nonequilibrium process in a finite-time schedule governed by the 
master equation. The left-hand side of JE is explicitly written as 

n-l 

(e- y ) ^ T = II {e-^^^^e^^l^^^Peq^ojio). (135) 

{a k } k=0 
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We take the first product of the above equation as, 

^|e-y(^;to) e 5tMo( < T 1 |ao;to)|p eq(ao;to) 

= P eq (<7i;ii)| 1 . (136) 

Repetition of the above manipulation in Eq. (135) yields the quantity in 
the right-hand side of JE as, 

n-l y y 

E^K^n^ 1 ^ ( 137 ) 

<y n fc=o 

where Z n = Zt- This is the case for a classical system on a heat bath, not 
for a quantum system. Although readers may think the above proof is not 
relevant to improvement of QAC, the formulation of JE for the classical 
system is available for QAC by aid of the classical-quantum mapping above 
introduced. 



4.2.3. Quantum Jarzynski annealing 

We here provide a novel protocol from the same spirit as JE by using 
the special quantum system (129). Let us consider to start from the trivial 
ground state with the uniform linear combination similarly to the case of the 
ordinary QA. This initial state expresses the high-temperature equilibrium 
state as |* eq (<o)) oc c^* ^ ^)/ 2 ^) with (3(t ) < 1. We introduce the ex- 
ponentiated pseudo work operator R(tk) = cxp(—Y(ak]tk)/'2). Observant 
readers might think it as a non-unitary operator, but we can construct this 
operation by considering an enlarged quantum system as detailed later. 
When we apply R(tk) to \^f eq (tk)) with the inverse temperature /3(tfc), the 
quantum state is changed into a state corresponding to the equilibrium dis- 
tribution with f3(tk+i). Then the application of the time-evolution operator 
U(cr'\<j;tk+i) = cxp(—i5tH q (a'\<j;tk+i)/h) does not alter the instantaneous 
quantum state, since it is the ground state of H q {a'\a\ tk+i). The resultant 
state after the repetition of the above procedure is 

n 

l*(*n)> « JJ{i2(tfc)^fc(<r fc |<r fc _i;tfc)}|*«i(to)>. 
fe=i 

(138) 

The product in the right-hand side is essentially of the same form as in 
Eq. (135). Instead of the exponentiated matrix of StM(ak+i\crk] tk), we use 
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the time-evolution operator Uk{<JkWk-i]tk) here. After the system reaches 
the state \^(t n )), we measure the obtained state by the projection onto 
a specified state a'. The probability is then given by |(c7'|^(t n ))| 2 , which 
means that the desired ground state can be obtained with the probability 
proportional to exp(— j3(t n )Ho), since \^(t n )) oc \^ eq (t n ))- If the above 
procedure continues up to /3(t n ) 1, the resultant wave function can yield 
the ground state of H - We call this procedure as the quantum Jarzynski 
annealing (QJA). 

We here emphasize the following three points. First, the protocol of QJA 
does not rely on the quantum adiabatic control. The computational time 
does not depend on the energy gap. In this sense, QJA does not suffer from 
the energy-gap closure differently from QAC. The required computational 
cost for realization of QJA is estimated from the number of the unitary 
gates as will be discussed below. Second, from the property of JE, the 
result is independent of the schedule to tune the parameter, T, in the above 
manipulations. Third, we do not need the repetition of the pre-determined 
process to deal with all fluctuations in the nonequilibrium-process average 
as in the ordinary JE, since the classical ensemble is mapped to the quantum 
wave function. In addition, if we obtain the final wave function, the output 
can give the ground state we desire with a very high probability since 
(3(t n ) » 1. a 

4.2.4. Problems in measurement of answer 

So far, so good. No problems seem to exist in the realization and perfor- 
mance of QJA. Unfortunately, we can find a serious problem to efficiently 
solve the hard optimization problem by QJA. In order to implement QJA, 
we must prepare a peculiar operator, the exponentiated pseudo work op- 
erator R(tk) = exp(— y(cr; tfc)/2), which looks like a non-unitary operator. 
We can realize this non-unitary operator for the original Hilbert space by 
preparing an enlarged quantum system with an ancilla qubit (another two- 
level quantum system) as I*, <j>\) = |>3/) £g> |^>i), where (pi is assumed to take 
and l. 44 We call the ancila qubit the computational state below. We ini- 
tially set 1^,01 =0). For simplicity, we restrict ourselves to the case with 



a Notice that several-time repetitions of experiments should be demanded since the out- 
put by quantum measurement is probabilistic. However we should emphasize that this 
point is not related with the theoretical property of JE attributed to rare events, neces- 
sity of all the realizations during the noncquilibrium process, but it comes from quantum 
nature. 
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Hq{o) > for any states. Let us introduce the following unitary operator 
for the enlarged quantum system as 

i^ fe ) = £k><<W V^^v/i-^MA 

= / CT ®r fc , (139) 

where j/fe(cr) = exp(— Yfc(cr; tfe)). We can obtain the weighted quantum sys- 
tem by applying this operator as y/yiJa)^, fa = 0). Then we regard 
R unlt (tk) as the exponentiated pseudo work operation R(tk) for the quan- 
tum state |*, fa = 0). The other probability amplitudes of i? unit |^, fa = 0) 
are given as 

(vM|i? unit -(t fc )|vM) = v ^M (140) 
(*,l|i? unit -(i fe )|vI/,0) = ^l-j/fcW. (141) 

The output in measurement of the quantum state, yT — ^(ct)^, fa = 1) 
is regarded as an undesired error state in our computation. Therefore we 
here do not employ such error states to find the ground state of the classical 
Hamiltonian H (a). Let us assume the uniform change of the inverse tem- 
perature (3(tk+i) — (3(tk) = Sfl for simplicity. In order to gain the relevant 
weight for obtaining the desired ground state of H , we must increase a pa- 
rameter corresponding to the inverse temperature up to (3(t n )e <~ 1, where 
e is the minimum energy gap of the "classical" Hamiltonian Ho (usually 
given by the energy unit). Therefore the number of steps of QJA is neces- 
sary up to n = (i(t n )/5(i <~ l/e8j3. The computational time n is the size 
of the enlarged quantum state we have to prepare as |\&, fa, ■ ■ ■ , (f> n ). The 
initial state of QJA in this case is l^eqC^o), 0, • • • ,0). The other states as 
1,0 ■ ■ ■ ,0), |*, 0, 1, • • • ,0), etc. such that several ancilla qubits take dif- 
ferent states from the initial ones as fa = 1 are regarded as the error states. 
To gain the weight up to exp(— f3(t n )H ), we perform the n-step exponen- 
tiated work operations as I„ <E> Y\ ® I2 <£> • • • <E> I n , I<j ® h ® Y2 ® ^3 <E> • • • ® I n , 
■ ■ ■ and I a ® I\ ® • • • ® 7„_i ® Y n , where Ij denotes the identity opera- 
tion. We can then obtain the desired state 0, 0, • • • ,0) with the weight 
as exp(—(3(t n )Ho). The weights for the other states, the error states, are 
given as (1 - exp(-/3(^) J ff )) 1/rl (exp(-/3(t„) J ffo)) 1 " 1/ " for 1,0, ■ ■ ■ ,0) 
and (l-exp(-/3(t„)i/ )) 2 /™(cxp(-/3(t rl )i/o)) 1 " 2/ " for |tf , 1, 1, 0, • • • , 0) and 
so on. 

Let us consider to obtain the meaningful outputs by the projective mea- 
surements. The desired state is simply one with fa = for any i. The 
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probability for obtaining the state \a, 0, 0, • • • ,0) is evaluated as 

e -P(t n )H (a) e -0(tn)H o (<r) 

= ELol 1 - e-' 3 ( t ") ff o(' T )) fc /"(e-' 9 (*")^( ff )) 1 - fe /" = 2^ ' 

(142) 

Thus the successful result is generated with the exponentially decreasing 
probability for increase of the system size N. We will demand more in- 
genious techniques to efficiently obtain the desired outputs. This is the 
remaining problem on QJA. 

4.3. Non-adiabatic quantum computation 

In the previous section, we consider to implement the manipulation in the 
left-hand side of JE in the quantum computation by use of the classical- 
quantum mapping. On the other hand, we have another direction of QA to 
improve the performance beyond QAC by directly considering the nonequi- 
librium behavior of quantum system. Again, in order to develop the theory 
of QA, we rely on the property of JE but for the quantum system. Instead of 
the adiabatic control, we consider to repeat non-adiabatic quantum anneal- 
ing (small or intermediate T) starting from a state chosen from equilibrium 
ensemble, not necessarily the ground state. We may not be able to easily 
reach the ground state of H by such processes, since the system does not 
keep the instantaneous ground state as in the adiabatic computation. We 
instead need to repeat the process many times to hit the ground state. In 
this way, the problem of long annealing time is expected to be replaced by 
many repetitions of non-adiabatic (possibly quick) evolution. We call such 
a procedure as non-adiabatic quantum annealing (NQA). 45 

4.3.1. Jarzynski equality for quantum system 

Before analysis on the detailed property of the non-adiabatic computation, 
we recall the Jarzynski equality but for isolated quantum system, 46 ' 47 while 
assuming its application to NQA. 

Let us consider to find the ground state as of the spin-glass Hamil- 
tonian H as in Eq. (127) by NQA. We prepare the dynamical quantum 
system following the time-dependent Hamiltonian (112). Initially we pick 
up a state from the canonical ensemble for H(0) = Hi = — ToEi^f an( ^ 
then let it evolve following the time-dependent Schrodingcr equation. The 
performed work in the isolated quantum system is given by the difference 
between the outputs of projective measurements of the initial and final 
energies, W — E m {T) — E n (0). Here m and n denote the indices of the 
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instantaneous eigenstates measured at the final and initial steps of NQA, 
H{T)\m{T)) = E m {T)\m(T)) and if(0)|n(0)) = £ n (0)|ra(0)), respectively. 
The time-evolution operator is given by the following unitary operator as 

U T = Tcxp jf dtH{t)^ , (143) 

where T denotes the time ordered product. Thus we can evaluate the tran- 
sition probability between the initial and final steps as 

Pm,n(0 -> T) = |(* ro (T)|f/ T |* n (0))| 2 . (144) 

Therefore the path probability for the noncquilibrium process starting from 
the equilibrium ensemble as 

WO^D ^W' , (145) 
z o(p; |%|) 

where we express the instantaneous partition function at each time t as 
Z t (f3; { Jij}). By directly evaluating the left-hand side of JE, we reach JE 
for the isolated quantum system as 



m.n ■> ' ' 



'VI 



m,n 

-PE m (T) 



(146) 



zomjij})' 

where we used the fact that the performed work W is a classical number 
and 

^P m , n (0-^T) = ^;<* m (T)|^r|*n(0)><*„(0)|4|* m (T)) 

n n 

= Y,{*m{T)\U T U^ m {T)) = 1. (147) 

m 

If we measure the physical observable Or at the last of the nonequilibrium 
process, we obtain another equation as 
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where the subscript on the square brackets in the right-hand side denotes 
the thermal average in the last equilibrium state with the inverse temper- 
ature /?. 

Below we show several observations by application of JE for the isolated 
quantum system to implementation of NQA. 

4.3.2. Performance of non-adiabatic quantum annealing 

First we discuss the possibility of NQA as a solver. Let us consider to 
measure equilibrium quantities through NQA. The ratio of Eqs. (146) and 
(148) gives 

(e-^) QA - (M9) 
The resultant equation suggests that the thermal average under the Hamil- 
tonian Hq on the right-hand side can be estimated by NQA on the left-hand 
side. This fact may be useful in the evaluation of equilibrium average, since 
the left-hand side is evaluated without slow adiabatic processes. In order 
to investigate the property of the ground state, we tune the inverse tem- 
perature into a very large value /3 > 1. We should be careful because the 
average on the left-hand side involves a non-extensive quantity, the expo- 
nentiated work, whose value fluctuates significantly from process to process. 
The average on the left-hand side must be calculated by many trials of an- 
nealing processes. Thus, rare events with large values of the exponentiated 
work (i.e. (i\W\ ^> To) would contribute to the average significantly, and 
we have to repeat the annealing process very many times in order to reach 
the correct value of the average. It usually needs very many, typically ex- 
ponentially many, repetitions. Thus the difficulty has not been relaxed yet 
in general, but the present new perspective may lead to different methods 
and tools than conventional ones to attack the problem. 

4.4. Analyses on non-adiabatic quantum annealing 

Unfortunately, we have not reached any answers on the performance of 
NQA. Instead We here evaluate several properties in nonequilibrium pro- 
cess as in NQA for the particular spin glasses. We can exactly analyze 
nonequilibrium behavior by combination of JE with the gauge transforma- 
tion, although, in general, there are few exact results in nonequilibrium 
quantum dynamical system with many components. 

Following the prescription of the Jarzynski equality, we consider a repe- 
tition of NQA starting from the equilibrium ensemble. Let us remember the 
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whole Hamiltonian of QA for the typical spin glasses. The initial Hamilto- 
nian is given only by the transverse field H(0) = H\, which means a trivial 
uniform distribution. Consequently, as a starting point of our analyses, we 
write down the specialized JE to the case for NQA as 



-fjW\ _ Z T(f3,{ J ij}) 

(2 cosh f3T o y 



(e-^)QA = Z^w - (150) 



We assume that the interactions {Jij} follow the distribution function for 
the ±J Ising model (29), which is better to be rewritten as 

where we do not use K = (3 J for transparency, and exp(— 2(3 P J) = (1 — p) /p. 



4.4.1. Gauge transformation for quantum spin systems 

For several special spin glasses as the ±J Ising model, the gauge trans- 
formation is available for analyses on the dynamical property even under 
quantum fluctuations. The time-dependent Hamiltonian as in Eq. (115) is 
invariant under the following local transformation, 

of -+ of, aV h. of -+ Jij -»• JiMi (Vi,j). (152) 

where £»(= ±1) is called as a gauge variable. This transformation is de- 
signed to preserve the commutation relations between different components 
of Pauli matrix. 48 We skillfully use the gauge transformation to analyze the 
dynamical property of the nonequilibrium behavior in NQA. 



4.4.2. Relationship between two different paths of NQA 

Below, we reveal several properties inherent in NQA by the gauge trans- 
formation. Let us take the configurational average of Eq. (150) over all the 
realizations of {Jij} for the special case with [3 — (i\ and f3 p — fj 2 as 



[<e-^) Q A] 



Zt(Pi;{J^}) 

(2coshf3 1 T ) N 



(153) 



The right-hand side is written explicitly as 



[<e >QA] ^ = .f . (2cosh/^ (2cosh/3 iro )- (154) 
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Let us here apply the gauge transformation and sum over all possible con- 
figurations of the gauge variables {&}■ We obtain, after dividing the result 
by2 N , 

r/ P -/w\ j _ z T(fa;{Jij})z T (fa;{Jij}) n __. 

L( >QaJ * " ^ 2^(2 cosh (2cosh/3 1 r ) JV - ( 55) 

l7ij J 

A similar quantity of the average of the exponentiated work for the spin 
glass with the inverse temperature fa and the parameter for the quenched 
randomness fa gives 

* " { X} (2cosh/3 2 r ) JV - ( 56) 

Comparing Eqs. (155) and (156), we find the following relation between 
two different non-adiabatic processes, 

17 -PiW\ i 17 1 /cosh/?iJ\ ^ B /cosh/3 2 r \ ^ 

We describe the two different paths of NQA related by this equality in 
Fig. 20. Setting fa = in Eq. (157), (implying p = 1/2, the symmetric 
distribution or the high-temperature limit), we find a simple equality on 
the performed work during NQA 

r/ -/W\ I (cosh ft J)** 

^ <e )qa J° = (cosher,,)"- (158) 

The symmetric distribution (/3 2 = on the left-hand side) makes it possible 
to reduce the right-hand side to the above trivial expression. It is remarkable 
that NQA, which involves very complex dynamics, satisfies such a simple 
identity irrespective of the speed of annealing T. If we apply the Jensen 
inequality to the above equality, we can obtain the lower bound for the 
performed work as 

Here we generalize the inverse temperature to j3 from the specific choice 
fa. This lower bound is loose, since the direct application of the Jensen 
inequality to JE for NQA yields, after the configurational average with the 
symmetric distribution, 

( N B \ 
'^f) ' <«»> 
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Fig. 20. Two different processes of NQA in Eq. (157). The left-hand side of Eq. 
(157) represents the annealing process ending at the upper-right black dot and 
the right-hand side terminates at the lower-left dot. Three phases expressed by 
the same symbols as in Fig. 12 are separated by solid curves and a vertical line. 
The dotted line expresses the NL /3 p = /3. 

where D(f3\j3') is the Kullback-Leibler divergence defined as 

D (m = £ P,{{j l3 })\o g p ^ {J ;^ (lei) 

Here we defined the marginal distribution for the specific configuration 
{ J^} summed over all the possible gauge transformations, 

Since the Kullback-Leibler divergence does not become non-negative, the 
work performed by the transverse field during a nonequilibrium process in 
in the symmetric distribution (i.e. the left-hand side of Eq. (160)) does not 
exceed the second quantity on the right-hand side of Eq. (160). This fact 
means Eq. (159) was looser. 
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4.4.3. Exact relations involving inverse statistics 

Beyond the above results, we can perform further non-trivial analyses for 
the nonequilibrium process in the special conditions. Let us next take the 
configurational average of the inverse of the Jarzynski equality, Eq. (150), 

as 



1 




"(2cosh/3r ) iV 


[(c- fjW ) Q A_ 


P 


_ Z T (P;{Jij}) 



(163) 



P P 



The application of the gauge transformation to the right-hand side yields 



1 



(e-^> 



QA 



E 



exp (/3 P E(y> JiMo) (2cosh/3r ) 



A" 



(2 cosh /3 P J) 



-Vi, 



(164) 



By summing the right-hand side over all the possible configurations of {&} 
and dividing the result by 2^, we reach 



1 



(e-^)QA 



E 



Z T {(3 p -{J l3 }) (2cosh/3r ) 



N 



0p 2^(2 cosh/V)^ ZtH3; {Jij}) ' 



If we set p p — p on the Nishimori line, this equation reduces to 



1 



(cosh /3To 



p (cosh/W)^' 



.<e-^) Q A. 

Comparison of Eqs. (158) and (166) leads us to 

1 



[<e-^>QA] 



(165) 



(166) 



(167) 



^.(e-^)QA. 

As shown in Fig. 21, two completely different processes are nontrivially 
related by the resultant relation: One toward the Nishimori line and the 
other for the symmetric distribution. 

Let us further consider the inverse of Eq. (148) for the two-point corre- 
lation Ot = a t aZ j- We take the configurational average of both sides under 
the condition f3 p = P as 



(afCTfe-^) QA 



N 



(cosh/3r ) 
(cosh PJ) n b 



(° z i° z i)p 



(168) 



The quantity on the right-hand side becomes unity by the gauge transfor- 
mation as has been shown in the literatures. 7 ' 22 We thus obtain a simple 
exact relation 

_ (cosh/JFo)^ 



(afaje-^) 



QA 



(cosh pj) 



N B 



(169) 
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•p 



Fig. 21. Two different nonequilibrium processes in NQA through Eq. (167). We 
use the same symbols as in Fig. 12. The white circle denotes the target of the 
process on the right-hand side of Eq. (167), whereas the black dot is for the 
left-hand side. 

which is another exact identity for processes of NQA. 

The importance of the above equalities is still not clear. However it is 
true that we have little exact results on the nonequilibrium behavior of the 
spin glasses driven by quantum fluctuations as the transverse field. When 
we realize the quantum spin systems in experiments, the above results can 
play a roll as their indicators to verify their precisions and conditions. We 
should emphasize that the specialized tool to analyze spin glasses has facili- 
tated the dynamical property in nonequilibrium process driven by quantum 
fluctuations. Beyond our analyses introduced here, we hope that more fas- 
cinating results on quantum computations would be obtained in the future 
with recourse to several techniques developed in statistical mechanics. 

5. Summary 

We looked over two topics lying between quantum information processing 
and statistical mechanics. 
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The first was the quantum error correction using the property of the 
topology. We prepare the redundant physical qubits in order to express the 
logical qubits we encode the specific information. Although we should deal 
with the quantum-many body systems, the used technique was proposed 
from the classical method in statistical mechanics. Statistical mechanics 
is available to facilitate to identify its precise location of the theoretical 
limitation to successfully infer the original state. In this direction with 
rich knowledge in statistical mechanics, we will be able to propose another 
ingenious way to give the quantum state more robustness and resilience as 
in the case of the depolarizing channel. 

The second part was to improve the performance of quantum annealing 
away from the adiabatic-control case. The theoretical key was the Jarzyn- 
ski equality. We suggest two ways to overcome the particular bottleneck 
in the adiabatic computation of quantum annealing. Both of the schemes 
would be needed for many-time repetition to produce the desired results. 
However we hope that the studies in this direction would give a novel tech- 
nique beyond the ordinary limitations. The key point will be to enhance 
the possibility with the desired states. In the classical counterpart, recently 
several techniques are proposed with use of the population to reach the 
desired distribution, that is equilibrium state. 49-52 Such a technique as in 
the classical case enables us to generate the particular quantum state with 
higher probability than expected by use of skillful techniques. 

The bridge between the quantum information processing and statistical 
mechanics continues to lead us to frontiers, where we encounter novel and 
surprising results beyond expected ones from knowledge we obtain only 
from each side. We are at a position to look at the birth of such a fascinating 
interdisciplinary field. Don't miss it. 
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